Genetic Mapping and Characterization of Verticillium Wilt Resistance in a Recombinant Inbred Population of Upland Cotton.

Verticillium wilt (VW) is an important and widespread disease of cotton and once established is long-lived and difficult to manage. In Australia, the non-defoliating pathotype of Verticillium dahliae is the most common, and extremely virulent. Breeding cotton varieties with increased VW resistance is the most economical and effective method of controlling this disease and is greatly aided by understanding the genetics of resistance. This study aimed to investigate VW resistance in 240 F7 recombinant inbred lines (RIL) derived from a cross between MCU-5, which has good resistance, and Siokra 1-4, which is susceptible. Using a controlled environment bioassay, we found that resistance based on plant survival or shoot biomass was complex but with major contributions from chromosomes D03 and D09, with genomic prediction analysis estimating a prediction accuracy of 0.73 based on survival scores compared to 0.36 for shoot biomass. Transcriptome analysis of MCU-5 and Siokra 1-4 roots uninfected or infected with V. dahliae revealed that the two cultivars displayed very different root transcriptomes and responded differently to V. dahliae infection. Ninety-nine differentially expressed genes were located in the two mapped resistance regions and so are potential candidates for further identifying the genes responsible for VW resistance.

Cotton Breeding in Australia: Meeting the Challenges of the 21st Century.

The Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program is the sole breeding effort for cotton in Australia, developing high performing cultivars for the local industry which is worth∼AU$3 billion per annum. The program is supported by Cotton Breeding Australia, a Joint Venture between CSIRO and the program's commercial partner, Cotton Seed Distributors Ltd. (CSD). While the Australian industry is the focus, CSIRO cultivars have global impact in North America, South America, and Europe. The program is unique compared with many other public and commercial breeding programs because it focuses on diverse and integrated research with commercial outcomes. It represents the full research pipeline, supporting extensive long-term fundamental molecular research; native and genetically modified (GM) trait development; germplasm enhancement focused on yield and fiber quality improvements; integration of third-party GM traits; all culminating in the release of new commercial cultivars. This review presents evidence of past breeding successes and outlines current breeding efforts, in the areas of yield and fiber quality improvement, as well as the development of germplasm that is resistant to pests, diseases and abiotic stressors. The success of the program is based on the development of superior germplasm largely through field phenotyping, together with strong commercial partnerships with CSD and Bayer CropScience. These relationships assist in having a shared focus and ensuring commercial impact is maintained, while also providing access to markets, traits, and technology. The historical successes, current foci and future requirements of the CSIRO cotton breeding program have been used to develop a framework designed to augment our breeding system for the future. This will focus on utilizing emerging technologies from the genome to phenome, as well as a panomics approach with data management and integration to develop, test and incorporate new technologies into a breeding program. In addition to streamlining the breeding pipeline for increased genetic gain, this technology will increase the speed of trait and marker identification for use in genome editing, genomic selection and molecular assisted breeding, ultimately producing novel germplasm that will meet the coming challenges of the 21st Century.

Genomic prediction of cotton fibre quality and yield traits using Bayesian regression methods.

Genomic selection or genomic prediction (GP) has increasingly become an important molecular breeding technology for crop improvement. GP aims to utilise genome-wide marker data to predict genomic breeding value for traits of economic importance. Though GP studies have been widely conducted in various crop species such as wheat and maize, its application in cotton, an essential renewable textile fibre crop, is still significantly underdeveloped. We aim to develop a new GP-based breeding system that can improve the efficiency of our cotton breeding program. This article presents a GP study on cotton fibre quality and yield traits using 1385 breeding lines from the Commonwealth Scientific and Industrial Research Organisation (CSIRO, Australia) cotton breeding program which were genotyped using a high-density SNP chip that generated 12,296 informative SNPs. The aim of this study was twofold: (1) to identify the models and data sources (i.e. genomic and pedigree) that produce the highest prediction accuracies; and (2) to assess the effectiveness of GP as a selection tool in the CSIRO cotton breeding program. The prediction analyses were conducted under various scenarios using different Bayesian predictive models. Results highlighted that the model combining genomic and pedigree information resulted in the best cross validated prediction accuracies: 0.76 for fibre length, 0.65 for fibre strength, and 0.64 for lint yield. Overall, this work represents the largest scale genomic selection studies based on cotton breeding trial data. Prediction accuracies reported in our study indicate the potential of GP as a breeding tool for cotton. The study highlighted the importance of incorporating pedigree and environmental factors in GP models to optimise the prediction performance.

Characterization and Genetic Mapping of Black Root Rot Resistance in Gossypium arboreum L.

Black root rot (BRR) is an economically important disease of cotton and other crops, especially in cooler regions with short growing seasons. Symptoms include black discoloration of the roots, reduced number of lateral roots and stunted or slow plant growth. The cultivated tetraploid Gossypium species are susceptible to BRR. Resistance to BRR was identified in G. arboreum accession BM13H and is associated with reduced and restricted hyphal growth and less sporulation. Transcriptome analysis indicates that BM13H responds to infection at early time points 2- and 3-days post-inoculation, but by day 5, few differentially expressed genes are observed between infected and uninfected roots. Inheritance of BM13H resistance to BRR was evaluated in an F6 recombinant inbred population and shows a single semi-dominant locus conferring resistance that was fine mapped to a region on chromosome 1, containing ten genes including five putative resistance-like genes.

Genetic mapping and transcriptomic characterization of a new fuzzless-tufted cottonseed mutant.

Fiber mutants are unique and valuable resources for understanding the genetic and molecular mechanisms controlling initiation and development of cotton fibers that are extremely elongated single epidermal cells protruding from the seed coat of cottonseeds. In this study, we reported a new fuzzless-tufted cotton mutant (Gossypium hirsutum) and showed that fuzzless-tufted near-isogenic lines (NILs) had similar agronomic traits and a higher ginning efficiency compared to their recurrent parents with normal fuzzy seeds. Genetic analysis revealed that the mutant phenotype is determined by a single incomplete dominant locus, designated N5. The mutation was fine mapped to an approximately 250-kb interval containing 33 annotated genes using a combination of bulked segregant sequencing, SNP chip genotyping, and fine mapping. Comparative transcriptomic analysis using 0-6 days post-anthesis (dpa) ovules from NILs segregating for the phenotypes of fuzzless-tufted (mutant) and normal fuzzy cottonseeds (wild-type) uncovered candidate genes responsible for the mutant phenotype. It also revealed that the flanking region of the N5 locus is enriched with differentially expressed genes (DEGs) between the mutant and wild-type. Several of those DEGs are members of the gene families with demonstrated roles in cell initiation and elongation, such as calcium-dependent protein kinase and expansin. The transcriptome landscape of the mutant was significantly reprogrammed in the 6 dpa ovules and, to a less extent, in the 0 dpa ovules, but not in the 2 and 4 dpa ovules. At both 0 and 6 dpa, the reprogrammed mutant transcriptome was mainly associated with cell wall modifications and transmembrane transportation, while transcription factor activity was significantly altered in the 6 dpa mutant ovules. These results imply a similar molecular basis for initiation of lint and fuzz fibers despite certain differences.

Genetic Identification and Transcriptome Analysis of Lintless and Fuzzless Traits in Gossypium arboreum L.

Cotton fibres, as single cells arising from the seed coat, can be classified as lint and fuzz according to their final length. Gossypium arboreum is a cultivated diploid cotton species and a potential donor of the A subgenome of the more widely grown tetraploid cottons. In this study, we performed genetic studies on one lintless and seven fuzzless G. arboreum accessions. Through association and genetic linkage analyses, a recessive locus on Chr06 containing GaHD-1 was found to be the likely gene underlying the lintless trait. GaHD-1 carried a mutation at a splicing acceptor site that resulted in alternative splicing and a deletion of 247 amino acid from the protein. The regions containing GaGIR1 and GaMYB25-like were found to be associated with fuzz development in G. arboreum, with the former being the major contributor. Comparative transcriptome analyses using 0-5 days post-anthesis (dpa) ovules from lintless, fuzzless, and normal fuzzy seed G. arboreum accessions revealed gene modules and hub genes potentially important for lint and fuzz initiation and development. Three significant modules and 26 hub genes associated with lint fibre initiation were detected by weighted gene co-expression network analysis. Similar analyses identified three vital modules and 10 hub genes to be associated with fuzz development. The findings in this study contribute to understanding the complex molecular mechanism(s) regulating fibre initiation and development and indicate that G. arboreum may have fibre developmental pathways different from tetraploid cotton. It also provides candidate genes for further investigation into modifying fibre development in G. arboreum.

Recapitulating early development of mouse musculoskeletal precursors of the paraxial mesoderm in vitro.

Body skeletal muscles derive from the paraxial mesoderm, which forms in the posterior region of the embryo. Using microarrays, we characterize novel mouse presomitic mesoderm (PSM) markers and show that, unlike the abrupt transcriptome reorganization of the PSM, neural tube differentiation is accompanied by progressive transcriptome changes. The early paraxial mesoderm differentiation stages can be efficiently recapitulated in vitro using mouse and human pluripotent stem cells. While Wnt activation alone can induce posterior PSM markers, acquisition of a committed PSM fate and efficient differentiation into anterior PSM Pax3+ identity further requires BMP inhibition to prevent progenitors from drifting to a lateral plate mesoderm fate. When transplanted into injured adult muscle, these precursors generated large numbers of immature muscle fibers. Furthermore, exposing these mouse PSM-like cells to a brief FGF inhibition step followed by culture in horse serum-containing medium allows efficient recapitulation of the myogenic program to generate myotubes and associated Pax7+ cells. This protocol results in improved in vitro differentiation and maturation of mouse muscle fibers over serum-free protocols and enables the study of myogenic cell fusion and satellite cell differentiation.

A Gradient of Glycolytic Activity Coordinates FGF and Wnt Signaling during Elongation of the Body Axis in Amniote Embryos.

Mammalian embryos transiently exhibit aerobic glycolysis (Warburg effect), a metabolic adaptation also observed in cancer cells. The role of this particular type of metabolism during vertebrate organogenesis is currently unknown. Here, we provide evidence for spatiotemporal regulation of glycolysis in the posterior region of mouse and chicken embryos. We show that a posterior glycolytic gradient is established in response to graded transcription of glycolytic enzymes downstream of fibroblast growth factor (FGF) signaling. We demonstrate that glycolysis controls posterior elongation of the embryonic axis by regulating cell motility in the presomitic mesoderm and by controlling specification of the paraxial mesoderm fate in the tail bud. Our results suggest that glycolysis in the tail bud coordinates Wnt and FGF signaling to promote elongation of the embryonic axis.

Global transcriptomic profiling in barramundi (Lates calcarifer) from rivers impacted by differing agricultural land uses.

Most catchments discharging into the Great Barrier Reef lagoon have elevated loads of suspended sediment, nutrients, and pesticides, including photosystem II inhibiting herbicides, associated with upstream agricultural land use. To investigate potential impacts of declining water quality on fish physiology, RNA sequencing (RNASeq) was used to characterize and compare the hepatic transcriptomes of barramundi (Lates calcarifer) captured from 2 of these tropical river catchments in Queensland, Australia. The Daintree and Tully Rivers differ in upstream land uses, as well as sediment, nutrient, and pesticide loads, with the area of agricultural land use and contaminant loads lower in the Daintree. In fish collected from the Tully River, transcripts involved in fatty acid metabolism, amino acid metabolism, and citrate cycling were also more abundant, suggesting elevated circulating cortisol concentrations, whereas transcripts involved in immune responses were less abundant. Fish from the Tully also had an increased abundance of transcripts associated with xenobiotic metabolism. Previous laboratory-based studies observed similar patterns in fish and amphibians exposed to the agricultural herbicide atrazine. If these transcriptomic patterns are manifested at the whole organism level, the differences in water quality between the 2 rivers may alter fish growth and fitness. Environ Toxicol Chem 2017;36:103-112. © 2016 SETAC.

Differentiation of pluripotent stem cells to muscle fiber to model Duchenne muscular dystrophy.

During embryonic development, skeletal muscles arise from somites, which derive from the presomitic mesoderm (PSM). Using PSM development as a guide, we establish conditions for the differentiation of monolayer cultures of mouse embryonic stem (ES) cells into PSM-like cells without the introduction of transgenes or cell sorting. We show that primary and secondary skeletal myogenesis can be recapitulated in vitro from the PSM-like cells, providing an efficient, serum-free protocol for the generation of striated, contractile fibers from mouse and human pluripotent cells. The mouse ES cells also differentiate into Pax7(+) cells with satellite cell characteristics, including the ability to form dystrophin(+) fibers when grafted into muscles of dystrophin-deficient mdx mice, a model of Duchenne muscular dystrophy (DMD). Fibers derived from ES cells of mdx mice exhibit an abnormal branched phenotype resembling that described in vivo, thus providing an attractive model to study the origin of the pathological defects associated with DMD.

A novel conserved mechanism for plant NLR protein pairs: the "integrated decoy" hypothesis.

Plant immunity is often triggered by the specific recognition of pathogen effectors by intracellular nucleotide-binding, leucine-rich repeat receptors (NLR). Plant NLRs contain an N-terminal signaling domain that is mostly represented by either a Toll-interleukin1 receptor (TIR) domain or a coiled coil (CC) domain. In many cases, single NLR proteins are sufficient for both effector recognition and signaling activation. However, many paired NLRs have now been identified where both proteins are required to confer resistance to pathogens. Recent detailed studies on the Arabidopsis thaliana TIR-NLR pair RRS1 and RPS4 and on the rice CC-NLR pair RGA4 and RGA5 have revealed for the first time how such protein pairs function together. In both cases, the paired partners interact physically to form a hetero-complex receptor in which each partner plays distinct roles in effector recognition or signaling activation, highlighting a conserved mode of action of NLR pairs across both monocotyledonous and dicotyledonous plants. We also describe an "integrated decoy" model for the function of these receptor complexes. In this model, a plant protein targeted by an effector has been duplicated and fused to one member of the NLR pair, where it acts as a bait to trigger defense signaling by the second NLR upon effector binding. This mechanism may be common to many other plant NLR pairs.

Next generation sequence analysis of the transcriptome of Sydney rock oysters (Saccostrea glomerata) exposed to a range of environmental stressors.

Sydney rock oysters (Saccostrea glomerata) were exposed to environmental stressors at contaminated field sites or in a controlled laboratory setting. RNA seq transcriptome data were generated for the gill and digestive gland using Roche's 454 pyrosequencing technology. 28,685 contigs were de novo assembled which encoded 11,671 different protein products. The data will act as a reference for future studies in ecology, immunology and environmental toxicology.

RNA-Seq analysis of the toxicant-induced transcriptome of the marine diatom, Ceratoneis closterium.

Diatoms are of enormous ecological importance as they account for as much as 20% of global primary production, yet they are still understudied from a genomic perspective. The benthic diatom Ceratoneis closterium is well-characterized from an ecotoxicological perspective including its use in ecotoxicological risk assessments and investigating the mode-of-action of metal toxicity. However, this organism has little sequence information available. In this study, 454 pyrosequencing of the stressor-responsive transcriptome was undertaken. These transcripts could be used to characterize general physiological processes such as photosynthesis and respiration, as well as to enable a description of the ecotoxicogenomic responses of this organism. After a 96 h exposure to the concentration of toxicant that inhibited growth rate by 10% (IC10) for the following common coastal contaminants: ammonia, copper, crude oil and simazine (a photosystem II inhibiting herbicide), diatom cells were harvested for RNA extraction and their transcriptomes characterized via 454 pyrosequencing. This resulted in 1.25 million reads, which were assembled into 4768 contigs, when contigs encoding rRNA were removed. More than 80% of the remaining contigs had an ortholog in the BLASTx protein databases. These contigs represented 1660 unique transcripts. The role of these transcripts in stress response, as well as photosynthesis and respiration is discussed. Overall, this study greatly enhances the genomic information available for this important taxonomic group.

454 pyrosequencing-based analysis of gene expression profiles in the amphipod Melita plumulosa: transcriptome assembly and toxicant induced changes.

Next generation sequencing using Roche's 454 pyrosequencing platform can be used to generate genomic information for non-model organisms, although there are bioinformatic challenges associated with these studies. These challenges are compounded by a lack of a standardized protocol to either assemble data or to evaluate the quality of a de novo transcriptome. This study presents an assembly of the control and toxicant responsive transcriptome of Melita plumulosa, an Australian amphipod commonly used in ecotoxicological studies. RNA was harvested from control amphipods, juvenile amphipods, and from amphipods exposed to either metal or diesel contaminated sediments. This RNA was used as the basis for a 454 based transcriptome sequencing effort. Sequencing generated 1.3 million reads from control, juvenile, metal-exposed and diesel-exposed amphipods. Different read filtering and assembly protocols were evaluated to generate an assembly that (i) had an optimal number of contigs; (ii) had long contigs; (iii) contained a suitable representation of conserved genes; and (iv) had long ortholog alignment lengths relative to the length of each contig. A final assembly, generated using fixed-length trimming based on the sequence quality scores, followed by assembly using the MIRA algorithm, produced the best results. The 26,625 contigs generated via this approach were annotated using Blast2GO, and the differential expression between treatments and control was determined by mapping with BWA followed by DESeq. Although the mapping generated low coverage, many differentially expressed contigs, including some with known developmental or toxicological function, were identified. This study demonstrated that 454 pyrosequencing is an effective means of generating reference transcriptome information for organisms, such as the amphipod M. plumulosa, that have no genomic information available in databases or in closely related sequenced species. It also demonstrated how optimization of read filtering protocols and assembly approaches changes the utility of results obtained from next generation sequencing studies, and establishes criteria to determine the quality of a de novo assembly in species lacking a reference genome. This new transcriptomic knowledge provides the genomic foundation for the creation of microarray and qPCR assays, serving as a reference transcriptome in future RNAseq studies, and allowing both the biology and ecotoxicology of this organism to be better understood. This approach will allow genomics-based methodology to be applied to a wider range of environmentally relevant species.